Scaling of a collapsed polymer globule in 2D 
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Extensive Monte Carlo data analysis gives clear evidence that collapsed linear polymers in two 
dimensions fall in the universality class of athermal, dense self-avoiding walks, as conjectured by 
B. Duplantier [Phys. Rev. Lett. 71, 4274 (1993)]. However, the boundary of the globule has self affirie 
roughness and does not determine the anticipated nonzero topological boundary contribution to 
entropic exponents. Scaling corrections are due to subleading contributions to the partition function 
corresponding to polymer configurations with one end located on the globule-solvent interface. 

PACS numbers: 05.70.Jk, 64.60.Ak, 36.20.Ey, 68.03.Cd 



Polymers in solution are the subject of intense studies 
since several decades Q. In the dilute case and in good 
solvent (high temperature T) excluded volume effects fa- 
vor swollen configurations for a long linear chain. On 
the other hand, in poor solvent (low T) effective attrac- 
tive interactions between monomers dominate and the 
typical conformations are those of a compact globule. 
The transition between swollen and collapsed polymer 
regimes is marked by the theta point 0, Q. A main 
achievement of the Coulomb gas Q and conformal in- 
variance approaches has been the exact characterization 
of the scaling properties of polymers in the swollen 0, [j| 
and theta regimes |a,Lll3] m two dimensions (2D). On the 
other hand, the situation is still far from settled as far as 
the entropic scaling of the collapsed phase is concerned. 

A crucial feature expected for a collapsed globule is the 
presence of a rather sharp boundary separating it from 
the surrounding solvent. This led Owczareck et al. to 
conjecture the presence of a factor growing like the expo- 
nential of the boundary length in the statistical partition 
sum descri bing the globule. It was then pointed out by 
Duplantier [lOj that with free boundary conditions such 
type of factor appears naturally for models of dense poly- 
mers (DP), in which a self avoiding walk (SAW) covers 
an assigned reg ion of the lattice visiting a fixed fraction 
of sites 0, Il2j . He also argued that a collapsed globule 
should have the exactly known scaling exponents of DP 
with smooth free boundary in 2D. Since DP have only 
excluded volume effects, this conjecture implies that col- 
lapsed configurations should not be sensibly influenced 
by the attractive interactions. A further, less obvious 
assumption is that the globule boundary acts as a sim- 
ple, smooth perimeter confining the polymer So far, 
the exponents conjectured in this way have never been 
confirmed by numerical investigations 0, 0] . 

In this Letter we show that the entropic scaling of col- 
lapsed polymers in 2D is consistent with the universality 
class of DP. Our analysis elucidates geometrical prop- 
erties of the globule-solvent interface and the role they 
play in determining exponents and strong finite size cor- 



rections to the asymptotic behaviors. 

Let us consider linear SAW's w of |to| = N steps 
on square lattice Q (Fig. QJ. If to each pair of near- 
est neighbor sites visited not consecutively by the SAW 
(contact) is associated an attractive potential energy — e 
(e > 0), the model undergoes the theta collapse transi- 
tion upon varying T. This transition can be monitored 
from the asymptotic behavior of the partition function 
Zn = Yj\w\=n exp[f where the sum extends to 

all possible w with an end at a fixed origin, and C(w) 
is the number of nearest neighbor contacts in w. Below 
the theta temperature T@ , Zn is expected [t| to have the 
following asymptotic behavior: 



Z N {T) ~Aii(T) N ni(T) vl "Ni 



(1) 



where A is an amplitude, /i(T) is a bulk free energy per 
step depending also on lattice structure, and ^i(T) < 1 
is a boundary term, associated with the existence of a 
sharp interface separating the typically globular region 
occupied by the SAW from the rest of the lattice. The 
factor i/n^ implies a boundary contribution to the di- 
mensionless total free energy, ln(Zjv). Indeed, y/N is the 
average number of SAW steps on the boundary of the 
globule, under the plausible assumption that this bound- 
ary has a fractal dimension equal to 1. If one restricts 
the sum in Eq. (1) to walks which start and end at the 
origin (polygons), the resulting partition (Zpj) has the 
same asymptotics as in Eq. (1), except for a different 
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FIG. 1: Collapsed SAW, with sites on the boundary (empty 
circles) distinguished from sites in the interior (dark circles). 
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exponent 70 replacing 7. Both 7 and 70 could take T- 
independent universal values in the collapsed polymer 
regime, as also implied by the conjecture in Ref. [Io| . 
On the basis of the analogy with DP [ill Il2| , Duplantier 
conjectured 7 — 70 = 19/16, and 70 = 5/6, implying 
7 = 97/46 10]. This value of 70 has a purely topologi- 
cal interpretation and was argued by assuming that for 
A — ► 00 the boundary of the globule becomes a smooth 
continuous arc without wedges [1(J| . An analysis of ex- 
tensive exact enumerations was performed in Ref. [lil ] , in 
order to determine 7 — 70 from extrapolations of Zjq jZ^ N . 
Considering such a ratio one expects the exponential and 
stretched exponential factors in Eq. (1) to simplify in 
numerator and denominator, leaving an A-dependence 
~ 7V 7 ~ 70 . The estimated 7 — 7n_~ 0.92 was inconsis- 
tent with the conjecture of Ref. [l0|. An effort to deter- 
mine 7 was subsequently made on the basis of extensive 
grand canonical Monte Carlo sampling : the estimate 
7 « 1.09, again appears to rule out the conjecture of 
Ref. [HI. 

A correct interpretation of this contradictory scenario, 
and a solution of the puzzle can be achieved once elu- 
cidated all physical consequences of the existence of a 
well defined interface between collapsed globule and sol- 
vent. This interface does not only imply the presence of 
a stretched exponential factor in the asymptotic parti- 
tion function in Eq. (1). Indeed, unlike in the swollen 
regime, the typical configurations of a polymer globule 
can be naturally partitioned into distinct groups, depend- 
ing on the location of the chain ends with respect to the 
boundary. This circumstance provides a natural source 
of scaling corrections to the asymptotic power law scal- 
ing in Eq. (1). For example, imagine one manages to re- 
strict the sum in Zn to compact configurations in which 
one end of the chain is located on the globule bound- 
ary, while the other one falls in the interior. We ask 
how the resulting "interior-boundary" partition function, 
Z$, should scale compared to Zn. The fact that scal- 
ing should not be influenced by attractive interactions 
and that the interfacial boundary has fractal dimension 
equal to 1, suggests that a simple A-dependent geomet- 
rical factor should connect the two partition functions: 
Z% ~ N~ X I 2 Z N . This factor A" 1 / 2 represents the prob- 
ability that one of the ends of a compact chain of A steps 
(globule area ~ A) ends on the globule boundary (length 
~ A 1 / 2 ), if one assumes that a chain end falls with equal 
probability anywhere within the globule. 

Of course, we also take into account that the fraction 
of walks with both ends on the boundary is in turn neg- 
ligible with respect to Zfy. Indeed, we can define also a 
restricted partition sum Z b N to which only chain config- 
urations with both ends on the boundary contribute. In 
this case the factor relating the restricted partition to Zjq 
will be A" 1 , since for large A one should regard as in- 
dependent events the occurrences of boundary locations 
for the two ends. Summarizing, if we further define Z l N 



as the partition restricted to configurations with both 
chain ends in the interior of the globule, the following 
asymptotics should be expected 

Z l N ~ A.^^N^ 1 (2a) 
Z% ~ Aib^p^N^- 1 (2b) 
Z b N ~ A b ^ N ^N^- 1 (2c) 

where Ai, A^, and A^ are suitable amplitudes, while 
"fib = 7 — 1/2 and 76 = 7 — 1. Clearly the sum of 
the three Z's above must yield Zn. So, the behaviors 
of Z$ and Z h Nl if confirmed, would also identify scaling 
corrections to the leading behavior in Eq. (1). The rel- 
ative magnitudes of the amplitudes will play a key role 
in determining up to what extent these corrections are 
important at finite N . 

In order to proceed we need a meaningful definition of 
the boundary and of the interior of a collapsed configura- 
tion. First, from now on we denote as neighbors the near- 
est and second neighbors of a site. For a SAW walk w, we 
define as boundary the set of visited sites which are neigh- 
bors of at least one non-visited site in communication 
with the exterior (Fig.^). In order to communicate with 
the exterior such non-visited site must be connected by 
at least one path of empty neighbor sites to the perimeter 
of a large lattice box enclosing the globule. The interior 
is then given by the sites visited by the walk which do 
not belong to the boundary. This boundary definition re- 
minds that of percolation cluster hull 0, H, ■ We were 
able to implement it by sampling long chain configura- 
tions via the nPERM algorithm (new pruned enriched 
Rosenbluth method) with importance sampling 0]. By 
this algorithm we evaluated weights proportional to Zn 
and to restricted partition functions defined above, up to 
7V max = 1920. The nPERM is an extremely efficient tool 
for sampling long compact SAW configurations and was 
even applied with success to study native structures of 
lattice proteins . The weight of the subset of configu- 
rations in which the SAW comes back to the origin gave 
an estimate of Z N . We explored the collapsed regime at 
e/T = 0.7, 0.77, and 0.85. For the case e/T = 0.85, on 
which we concentrated most efforts, we sampled 2 x 10 9 
linear chain configurations for each A (2 x 10 6 completely 
independent). This took more than one year CPU time 
on 2 GHz machines. 

We first checked how the average number, B^, of sites 
on the boundary [Fig. IHa)] grows with A. In Fig. G[a) 
we show Bm vs N in a log-log scale for e/T = 0.85. The 
data give clear evidence that asymptotically Bjq ~ A 1 / 2 . 
Thus, it makes sense to check if the scalings in Eqs. © 
are consistent with the data, and also provide the lead- 
ing scaling correction mechanism to Eq. (1). We consid- 
ered the ratios of all partition functions with respect to 
Z N . These ratios are reported as a function of A in the 
log-log plots of Fig. |3 which give an instructive repre- 
sentation of the role played by the scaling corrections at 
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FIG. 2: (a) Bn vs N for e/T = 0.85: data approaches a scal- 
ing ~ N 1 ^ 2 (dashed line) for large N. (b) Effective exponent, 

i.e. the slope of the data in (a), as a function of iV . The 
line is a linear fit extrapolating to ~ 0.5 for N —* oo. 



various N. The logarithmic slope of Zn/Z n vs N, which 
should be 7—70, appears to have strong corrections which 
still sensibly bend the curve at the highest values of N. 
Remarkably, the slope of Z % N /Z N approaches nearly the 
same value in this extremal region, suggesting that in- 
deed the leading source of corrections is primarily in Zj§, 
as discussed above. The fact that up to N « 1300 the 
curve for Z N remains below that for 27jy is due to a large 
ratio An, jAi . The extrapolated slopes of the curves are 
consistent with 7 — 7^ = Jib — 76 = 1/2 according to 
Eqs. © (see Fig.0J}. 

The above results suggest that the asymptotics of the 
collapsed regime should be best determined by isolating 
and studying Z l N . To confirm this, we evaluated an ef- 
fective entropic exponent 7 — 70 using a weighted linear 
least square fit of log 10 (ZN /Z N ) vs log 10 N (and similarly 
for Z l N ). To quantify the approach to the asymptotic 
scaling, we consider subsets of four consecutive points 
(Ni, N 2 , N 3 , Ni) and fit their slope in Fig. El The re- 

suits are plotted as a function of 1/N in Fig.0] where 
N = (NiN 2 N 3 N 4: ) 1 ^ 4 . In the case of Z N there is much 
curvature in the plot. Quite remarkably in the case of 
Z % N the plot appears pretty linear and extrapolates very 
close to 7 — 70 = 19/16 for N — > 00. This linearity 
clearly indicates that A r_1//2 represents the leading scal- 
ing correction. One can push further the analogy with 
DP, imagining that the N' = N — Bn monomers in the 
interior of the globule are representative of the DP phase, 
the remaining Bn ones just forming a kind of box bound- 
ary. In other words, the partition function can also be 
parametrized by the number N' of monomers in the in- 
terior, i.e. Zn> = Zpf. Repeating the same analysis 
with N' in place of N, we obtain effective exponents for 
Z N ,/Zjy, that lie on an almost horizontal straight line 
(Fig. 0} , showing that the scaling correction amplitude 
has been drastically reduced. This allows a more accu- 
rate extrapolation of 7 — 70 = 1.20(3), very close to the 
value 19/16 expected for DP. 
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FIG. 3: Partition function ratios for e/T = 0.85. For N < 100 
the dominant contribution to Zn comes from Z%, for 100 < 
N < 1300 Z* dominates, while Z l N becomes the leading term 
only for N > 1300. Lines emphasize the trends. 
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FIG. 4: Effective 7 — 70 as a function of N (dotted) or 

W~ 1/2 (dashed) for e/T = 0.85. Fits of 7^-70 and of ~/ b -~f 
are also shown. Arrows mark the DP value 7 — 70 = 19/16, 
and 7 — 70 — 1/2 and 7 — 70 — 1. 



As a further step we determined individual 7's by di- 
rectly fitting In Zn with a function of the form 



N In n + VN In fn + (7 - 1) In N + la A + A ± /VN (3) 

and Z N with a similar one [13. The term A 1 /^N corre- 
sponds to the expected scaling correction. With or with- 
out this term, we found that upon removal of the data at 
smaller ATs the fits are not always stable, again a signal 
of the presence of a strong correction to scaling. On the 
other hand, with the Ai/^/N term, the best asymptotic 
estimate of 7 is expected when all data are included. On 
this basis we estimated 7 = 1.18(4) and 70 = 0.02(7), 
again for e/T = 0.85. We also fitted data for e/T = 0.70. 
At this temperature one observes theta point values of 
7 and 7q. This means that this temperature is still too 
close to the critical value e/Te = 0.665 H3 to observe 
the DP scaling in chains with N < 2000. However, for 
e/T = 0.77 one finds results pretty consistent with those 
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TABLE I: Estimated and exact exponents. 



e/T 




walk: 7 




polygon: 70 




estim. 


exact 




estim. exact 


0.7 


1.11(3) 


8/7 " ~ 


1.14 


-0.15(5) -1/7 a ~ -0.14 


0.77 


1.20(5) 


19/16 b 


~ 1.19 


0.00(5) b 


0.85 


1.18(4) 


19/16 6 


~ 1.19 


0.02(7) b 



"Exact entropic exponent s of the theta point in 2D [3. 
'Exponents of DP in 2D llll . assuming 70 = 0. 



found for e/T = 0.85, confirming the expectation that 
the collapsed phase is described by T-independent ex- 
ponents. Our determinations are reported in Table I. 
Previous 7 estimates without scaling corrections in- 
cluded 7 = 1.09(8) at e/T ~ 0.765 and 7 = 1.11(6) at 
e/T ~ 0.788. While not testing as deeply the collapsed 
phase, these less sharp determinations are still compati- 
ble with ours. 

While the direct fits fully confirm the difference 7 — 70 
extrapolated above, consistent with the DP value 19/16, 
70 = 5/6, as conjectured in Ref. 0, appears definitely 
excluded. The fits suggest that 70 should be a much 
smaller number, possibly zero. The value 70 = 5/6 
was predicted in Ref. [Trj by assuming that the globule- 
solvent interface can be assimilated to a smooth wall 
without wedges in the DP model. However, the circum- 
stance that the boundary has a fractal dimension equal 
to 1, as directly verified here, does not rule out other pos- 
sibilities. As we argue below, one should expect a rough, 
sclf-affine boundary. The attractive forces guaranteeing 
the cohesion of the globule should determine a line ten- 
sion (~ ln/ii) for the boundary, which is probably the 
most important factor controlling its length. The fluc- 
tuations of the length of the a boundary should have a 
self-affine geometry, so that its average width grows like 
S^r ~ iV-/ 2 , C < 1 being the roughness exponent [l9|. In 
2D the roughness exponent of a line tension controlled 
boundary has £ = 1/2. One can argue the roughness ex- 
ponent C of the globule boundary again on the basis of 
the leading scaling correction ~ N^ 1 / 2 identified in the 
problem. If the boundary is sclf-affine, Bn should have a 
subleading correction ~ AT^ -1 , with positive amplitude, 
due to the rate of growth with A~ of the average width of 
the boundary profile. Thus, in our case we should expect 
£ = 1/2. This correction is clearly shown by our plot 

in Fig. HJb), where the effective exponents of the scal- 

1/2 

ing of B]s are plotted as a function of N : they vary 
1/2 

linearly as a function of N , with positive slope, and 
extrapolate to 0.507(3) as # -> 00 (dot-dashed line in 
the figure). A self-affine curve is not differentiable and it 
appears unjustified to represent the globule boundary as 
equivalent to a smooth contour for a DP 0] . The theory 
of DP [HI flij l , while still valid for the collapsed globule, 
is not applicable in the form appropriate for polymers 



within smooth boxes. Thus the topological argument 
leading to 70 = 5 /6 for a collapsed globule does not hold. 
One should expect a smaller value of 70, as the theory 
predicts that 70 is maximal for a DP with smooth bound- 
ary ^(J. Indeed, here we find 70 ~ 0, possibly coinciding 
exactly with 0, compatible with the hypothesis that a 
collapsed polymer is statistically equivalent to a DP with 
self-affine rough boundary. 

In summary we gave for the first time solid evidence 
that the entropic scaling of a collapsed polymer globule 
in 2D falls in the universality class of athermal DP. The 
existence of the solvent-globule interface appears crucial 
in several respects. Besides giving rise to the strong scal- 
ing corrections which hindered so far the analysis of the 
problem, with its nontrivial, self-affine stochastic geome- 
try the interface also determines an unexpected, close to 
zero value of the 70 exponent of collapsed rings. 
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